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Abstract 

Difference arrays are used in applications such as software testing, authentication codes and data compres¬ 
sion. Pseudo-orthogonal Latin squares are used in experimental designs. A special class of pseudo-orthogonal 
Latin squares are the mutually nearly orthogonal Latin squares (MNOLS) first discussed in 2002, with gen¬ 
eral constructions given in 2007. In this paper we develop row complete MNOLS from difference covering 
arrays. We will use this connection to settle the spectrum question for sets of 3 mutually pseudo-orthogonal 
Latin squares of even order, for all but the order 146. 
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1 Introduction 


Difference matrices are a fundamental tool used in the construction of combinatorial objects, generating a 
significant body of research that has identified a number of existence constraints. These difference matrices 
have been used for diverse applications, for instance, in the construction of authentication codes without secrecy 
HU, software testing and data compression [lOl . This diversity of applications, coupled with existence 

constraints, has motivated authors to generalise the definition to holey difference matrices, difference covering 
arrays and difference packing arrays, to mention just a few. 

In the current paper we are interested in constructing subclasses of cyclic difference covering arrays and 
exploiting these structures to emphasize new connections with other combinatorial objects, such as pseudo- 
orthogonal Latin squares. We use this connection to settle the existence spectrum for sets of 3 mutually 

pseudo-orthogonal Latin squares of even order, in all but one case. We begin with the formal definitions. 

A difference matrix (DM) over an abelian group (G, +) of order n is defined to be an nxk matrix Q = [q(i, j)] 
with entries from G such that, for all pairs of columns 0 < j, j' < k — 1, j ^ j ', the difference set 

A jj' = {q(i,j) - q{i,j') | 0 < i < n - 1} 

contains every element of G equally often, say A times. (See, for instance, [6], [S] and [9].) Note that we label the 
rows from 0 to n — 1 and the columns 0 to k— 1. Also to be consistent with later sections involving Latin squares 
and covering arrays our definition uses the transpose of the matrix given in [S] and (Sj. Since the addition of a 
constant vector, over G, to all rows and a constant vector to any column does not alter the set A j j >, we may 
assume that one row and one column contain only 0, the identity element of G. More precisely, to simplify later 
calculations, we will assume that all entries in the last row and last column of Q are 0. A difference matrix will 
be denoted DM(n, k ; A). If (G, +) is the cyclic group we refer to a cyclic difference matrix. 

Theorem 1.1. /(J Thm 17.5, p 411] A DM(n,k; A) does not exist if k > An. 

In the main, we will use difference matrices with k = 4, A = 1 and where possible we will work with cyclic 
difference matrices. In Section [3] we list a number of existence results that will be relevant to the current paper. 

A holey difference matrix (HDM) over an abelian group (G, +) of order n with a subgroup H of order h 
is defined to be an (n — h) x k matrix Q = [q(i, j)] with entries from G such that, for all pairs of columns 
0 < j,j' < k — 1, j ^ j ', the difference set 

A jj> = {q{i,j) ~ q{i,j') \ 0<i<n — h — l} 

contains every element of G \ H equally often, say A times. A holey difference matrix will be denoted 
HDM(fc,n;ft), where |G| = n and \H\ = h. If G is the cyclic group then we refer to a cyclic holey differ¬ 
ence matrix. 

Remark 1.2. As before a constant vector may be added to any column without affecting A jj> so we may assume 
that all entries in the last column of Q are equal to 0. However since H is a subgroup, 0 belongs to the hole. 
Consequently 0 does not occur in A jj' , and thus there will be no row containing two or more 0’s. Further since 
Aj k-i = G \ H, 0 < j < k — 2, the entries of H do not occur in the first k — 1 columns of Q. 

A difference covering (packing) array over an abelian group (G, +) of order n is defined to be an 77 x k matrix 
Q = [q(i, j)} with entries from G such that, for all pairs of distinct columns 0 < j, j' < k — 1, the difference set 

A j tj > = { q{i,j ) - \ 0 < i < g - 1} 

contains every element of G at least (at most) once. (See, for instance, [19] and [20].) A difference covering 
array will be denoted DCA(fc, ry; n) and a difference packing array will be denoted DPA(fc, 77 ; n). If (G, +) is the 
cyclic group, then the difference covering (packing) array is said to be cyclic. 

Difference covering arrays have been studied in their own right and are related to mutually orthogonal partial 
Latin squares and transversal coverings, with applications in information technology, see [T] and US¬ 
As before we may assume that the last row and last column of a DCA(fc, 77 ; n) contain only 0. 

In the papers m and [20], Yin constructs cyclic DCA(4, n + 1; n) for all even integers n, with similar results 
for cyclic difference packing arrays. Yin documents a number of product constructions for difference covering 
arrays, some of which will be reviewed in Section [3] and then adapted to construct difference covering arrays 
with specific properties; properties that build connections with pseudo-orthogonal Latin squares. 
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The additional properties that we seek are that 0 (the entry relating to identity element of G) occurs at 
least twice in each column of the DCA (k,n + l;n) and for pairs of columns, not including the last column, 
the repeated difference is not the element 0. Formally we are interested in DCA(fc,n + 1 ;n), Q = [q(i,j)], 
(0 < * < n, 0 < j < k — 1) satisfying the properties: 

PI. the entry 0 £ G occurs at least twice in each column of Q , and 

P2. for all pairs of distinct columns j and j', j ^ k — 1 ^ j', Ajjt = {q(i , j) — q(i, j') | 0 < i < n — 1} = G\ {0}, 

Note that this last property implies that Ajji contains a repeated difference that is not 0. 

The following example, of cyclic DCA(4, 7; 6) that satisfies PI and P2, is taken from (13]. 


B t 


0123450 

1350240 

3041520 

0000000 


In the next lemma we show that if G is the cyclic group over Z ra , then these conditions imply that for all 
distinct columns j and j', j ^ k — 1 ^ j’, 


A jt j> = {0,1,2,..., n/ 2, n/2,..., n - 1} 


with repetition retained. 

Lemma 1.3. If there exists a cyclic DCA{k,n + 1 ;n), Q = [q{i, j)\, (0 < i < n, 0 < j < k — 1) satisfying 
Properties PI and P2, then n is even. Further, given do such that do = q(i,j) — q{i,f) = q(i',j) ~ q(i',j'), for 

1 ^ i’ and k — 1 ^ j ^ j' ^ k — 1, then do = n/2. 

Proof. Let Q = [q(i,j)\ (0 < i < n, 0 < j < k — 1) represent the difference covering array. The definition 
requires that Z ra C A h j> and since column k — 1 of Q contains all zeros, Property PI implies that the remaining 
columns are permutations of the multi-set {0, 0, 1 , 2 ,..., n — 1}. 

Let do £ Z n \{0} represent the repeated difference in A jj>. Suppose n is odd and, without loss of generality, 
that column 0 is in standard form. Then, for all 0 < j < k — 2, ^17=0 Q(hl) = and 

v-4,. ,. ... (n — l)n 

= —2 - hd ° modn - 

i=0 

Consequently 2do = n(2u — n + 1), or equivalently n|2d 0 . But since n is odd, this leads to the contradiction, 
do £ Z n and n\do- Thus n is 2 p for some integer p , where p divides do, implying do = p. 

□ 

The remainder of this paper is organised as follows. In Section [2] we will draw the connection between 
DCA(fc, n + 1; n) and sets of mutually pseudo-orthogonal Latin squares and for a subclass of squares settle the 
spectrum question for all but a single order, namely 146. In Section[3]we review some of the general constructions 
for difference covering arrays and show that these constructions can be used to construct DCA(fc,n + 1 ;n) 
that satisfy Properties PI and P2. In Section [I] we give three new constructions for DCA(4, n + l;n)’s and 
consequently new families of mutually pseudo-orthogonal Latin squares. 

The notation [a, b ] = {a, a + 1,..., 6 — 1, 6} refers to the closed interval of integers from a to b. 

2 Pseudo-orthogonal Latin squares and difference covering arrays 

In this section we verify that cyclic difference covering arrays can be used to construct pseudo-orthogonal Latin 
squares. 

A Latin square of order n is an n x n array in which each of the symbols of Z n occurs once in every row and 
once in every column. Two Latin squares A = [a(i,j)\ and B = [b(i, j)], of order n , are said to be orthogonal if 

O = {{a(i,j), b(i,j)) | 0 < i, j < n - 1} = Z n x Z„. 
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A set of t Latin squares is said to be mutually orthogonal, f-MOLS(n), if they are pairwise orthogonal. A set of 
t idempotent MOLS(n), denoted f-IMOLS(n), is a set of t-MOLS(n) each of which is idempotent; that is, the 
cell (i, i) contains the entry i, for all 0 < i < n — 1. 

It is well known that difference matrices can be used to construct sets of mutually orthogonal Latin squares, 
see for instance j9j Lemma 6.12]. 

While the applications of orthogonal Latin squares are well documented, there are still many significant 
existence questions unanswered. For instance, it is known that there is no pair of MOLS(6), however it is 
not known if there exists a set of three MOLS(IO), or four MOLS(22), see [6]. The existence of a set of four 
MOLS(14) was established by Todorov [16] in 2012, but it is not known if there exists a set of five MOLS(14). 
Many of the existence results have been obtained using quasi-difference matrices or difference matrices with 
holes, see [6;. 

The importance and applicability of MOLSs combined with these difficult open questions has motivated 
authors, such as Raghavarao, Shrikhande and Shrikhande m and Bate and Boxall [2], to slightly vary the 
orthogonality condition to that of pseudo-orthogonal. A pair of Latin squares, A = [a(i,j)\ and B = [b(i,j)], of 
order n, is said to be pseudo-orthogonal if given O = {(a(i,j), b(i,j)) | 0 < i, j < n — 1}, for all o£Z„ 

I {(a,b(i,j)) | ( a,b(i,j )) € 0}\ = n- 1. 

That is, each symbol in A is paired with every symbol in B precisely once, except for one symbol with which it 
is paired twice and one symbol with which it is not paired at all. A set of t Latin squares, of order n , are said 
to be mutually pseudo-orthogonal if they are pairwise pseudo-orthogonal. 

The value and applicability of pseudo-orthogonal Latin squares has been established through applications 
to multi-factor crossover designs in animal husbandry p] , and strongly regular graphs [3] (though the definition 
varies here). Mutually nearly orthogonal Latin squares (MNOLS) are a special class of pseudo-orthogonal Latin 
squares, in that the set O does not contain the pair (a, a), for any a£Z„. Mutually nearly orthogonal Latin 
squares (MNOLS) were first discussed in a paper by Raghavarao, Shrikhande and Shirkhande in 2002 [15]. 

A natural question to ask is: Can we use difference techniques to construct mutually pseudo-orthogonal 
Latin squares? Raghavarao, Shrikhande and Shirkhande did precisely this and constructed mutually pseudo- 
orthogonal Latin squares from cyclic DCA(fc,n + 1 ;n) termed (k, n)-difference sets in [13]. The Raghavarao, 
Shrikhande and Shirkhande result is as follows. 

Theorem 2 . 1 . If there exists a cyclic DCA(t + 1, 2p + 1; 2p), Q' = that satisfies PI and P2, then there 

exists a set oft pseudo-orthogonal Latin squares of order 2 p. 

Proof. Recall that without loss of generality we may assume that the last row and column of Q' contain all 
zeros. Construct a new matrix Q = [q(i, j)] by removing the last row and last column from Q' and define a set 
of t arrays, L s = [l s (i,j)], 0 < s < t — 1, of order 2 p, by 

h(i,j) = q(i,s) +j(mod 2p), 0<i,j <2p-l. (2.1) 

It is easy to see that each column of L s is a permutation of Z 2 p and so L s is a Latin square. By Lemma 11.31 

A o,j’ = Wihj) - I 1 < * < 2 p} = (Z 2 p \ {0}) U {p} 

implying that when any two Latin squares are superimposed we obtain the set of ordered pairs ({Z 2 P x Z 2 P } \ 
{(a;, x) | 0 < x < 2p — 1}) U {(x, x + p) | 0 < x < 2p — 1} with repetition retained. 

□ 

If there exists a pair of pseudo-orthogonal Latin squares generated from cyclic difference covering arrays 
satisfying PI and P2, then there exists a pair of nearly orthogonal Latin squares. Conversely, a pair of nearly 
orthogonal Latin squares are necessarily pseudo-orthogonal Latin squares. Given this and the strong connection 
with papers m and [13] we will state all results in terms of mutually nearly orthogonal Latin squares. 

Raghavarao, Shrikhande and Shirkhande established bounds on the maximum number of Latin squares in 
a set of mutually nearly orthogonal Latin squares. This result provides bounds on k for DCA(fc, n + 1; n) that 
satisfy PI and P2. 

Lemma 2.2. Let p >2 be a positive integer. If there exists a cyclic DCA(k + l,2p + l;2p) that satisfies PI 
and P2, then k < p + 1. Further if p is even and there exists a DCA(k , 2p + 1; 2 p), then k < p + 1. 
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Proof. If k > p + 1 then there exists a set of more than (p + l)-MNOLS(2p), which contradicts Raghavarao, 
Shrikhande and Shrikhande result. Similarly for the second statement. 

□ 


2.1 Some interesting facts 

It is also interesting to note that the MNOLS(2p) constructed from cyclic difference covering arrays are essentially 
copies of the cyclic group. Consequently, these Latin squares are all bachelor squares , in that they have no 
orthogonal mate. 

We also note that these sets of Latin squares are row complete. A row complete Latin square, L = [, l(i,j )] 
is one in which the columns can be reordered in such a way that the set {(/(*, j + 1) | * £ Z ra ,0 < j < 
n — 1} = Z„ x Z n . So the set of entries obtained by taking pairs of adjacent cells in the same row, for all rows, 
gives the set of all ordered pairs on Z„. 

Williams [TS] verified that the columns of the Latin square corresponding to the Cayley table of the cyclic 
group can be rearranged to obtain a row complete Latin square. 

Each of the k MNOLS(2p) constructed from cyclic DCA(fc + 1, 2p + 1; 2 p) can be obtained by reordering the 
rows of the Cayley table of the cyclic group, without touching the columns. Hence simultaneously reordering 
the columns of these nearly orthogonal Latin squares will also produce row complete pseudo-orthogonal Latin 
squares. 


3 The spectrum for sets of 3 mutually nearly orthogonal Latin 
squares 

In 2007, Li and van Rees m continued the study of 3-MNOLS(n) and conjectured that they exist for all even 
n > 6. In a partial answer to this question, Li and van Rees proved the existence for small orders and orders 
greater than 356, (see also Fl2l). 

Theorem 3.1. \11\ Thm f.8] If 2p > 358, then there exists a 3-MNOLS(2p). 


This work was extended in 2014, when Demirkale, Donovan and Khodkar [7j developed further constructions 
for cyclic DCA(4, 2p + 1; 2 p) proving: 

Theorem 3.2. J?]/ There exist 3-MNOLS(2p), where 2p = 14, 22, 38, 46 mod 48. 

The next result lists known values for cyclic DCA(4,2 p + 1; 2 p) satisfying PI and P2, with 2 p < 356. 

Lemma 3.3. There exists cyclic DCA( 4, 2 p + 1; 2 p) for 2p = 6,8,..., 20, 22, 38, 46, 62, 70, 86, 94, 110, 118, 
134, 142, 158, 166, 182, 190, 206, 214, 230, 238, 254, 262, 278, 286, 302, 310, 326, 334, 350. 

Proof. The existence of orders 6 and 8 was given in {13; and orders 10, 12, 14, 16, 18, 20 in m- All the 
remaining cases were shown to exist in [7]. 

□ 

van Rees recently summarised these results and indicated that 3-MNOLS(2p) exist for all orders except 
possibly those given below. 

Lemma 3.4. [T7 | / A set of 3-MNOLS(2p) exists except possibly when 2p = 24, 26, 28, 30, 34, 36, 42, 50, 52, 54, 
58, 66, 74, 82, 92, 102, 106, 114, 116, 122, 124, 130, 138, 146, 148, 170, 172, 174, 178. 

In this section we will settle the existence question for all even orders except 2 p = 146. For completeness 
we list all values less than 2 p = 358 and, given the connection with row complete Latin squares, where possible 
we will use cyclic difference covering arrays to construct the Latin squares. Therefore we begin this section by 
reviewing relevant results from [8], [19] and [20] , and adapting these to construct cyclic DCA(4, 2?n + l;2m) 
that satisfy PI and P2. Where appropriate we will indicate how these results can be used to settle the spectrum 
question. 

We begin with the following straight forward result that is analogous to [8] Lem 2.3]. 
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Lemma 3.5. Suppose that there exists a HDM{k,n\h) over the group (G,+) with a hole over the subgroup H. 
Further suppose there exists a DCA(k , h+ 1; h) over H satisfying PI and P2. Then there exists a DCA{k , n+ 1; n ) 
over G satisfying PI and P2. Further suppose that the HDM(k,n\h) and DCA(k, h + l;h) are cyclic. Then 
there exists a cyclic DCA(k,n + l;n) satisfying PI and P2. 

Proof. In the cyclic case, let A = [a(i, j)] (0 < i < n — 1 — h, 0 < j < k — 1) represent the cyclic HDM(fc, n; h) 
and B = [b(i,j)\ (0 < i < h, 0 < j < k — 1) represent the cyclic DCA (k,h + 1 ;h). The definition of cyclic 
implies that H = {0, u, 2 u,.. .(h — l)n}, where n = uh and the proof of Lemma 11.31 implies that h is even and 
the repeated difference in , j k — 1 ^ j', of B is hu/2 = n/2. 

Set Q = [q(i,j)\ (0 < i < n, 0 < j < k) to be the concatenation of A with B and we obtain a cyclic 
DCA(fc,n + l;n) that satisfies PI and P2. 

The non-cyclic case follows similarly. 

□ 

Next we give a general product type construction taken from [8] and adapt it to construct cyclic difference 
covering arrays that satisfy PI and P2. 

Lemma 3.6. /3J Lem 2.6] If both a cyclic HDM(k,n; h) and a cyclic DM(n',k\ 1) exist, then so does a cyclic 
HDM(k,nn';hn'). In particular, if there exists a cyclic DM(n,k; 1) and a cyclic DM(n',k\ 1) then there exists 
cyclic DM(nn !, fc; 1). 

The first statement of Lemma T3.6I coupled with Lemma T3. 5 1 leads to the following straightforward result. 

Corollary 3.7. Suppose that there exists a cyclic HDM(k , n ; h), a cyclic DM(n ', k ; 1) and a cyclic DCA(k , hn'+ 
1 ;hn') that satisfies PI and P2. Then there exists a cyclic DCA(k,nn' + 1 ;nn') that satisfies PI and P2. 

The second statement of Lemma 13.61 can also be adapted. 

Lemma 3.8. Suppose a cyclic DM(n,k\T), a cyclic DM(n',k; 1) and a cyclic DCA{k,n' + l;n') satisfying PI 
and P2 exist. Then there exists a cyclic DCA{k , nn' + 1; nn') that satisfies PI and P2. 

Proof. This result can be obtained by taking a hole of size 1 in Corollary 13.71 or as follows. Let A = [a(i,j)] 
(0 < i < n— 1, 0 < j < k— 1) represent the cyclic DM(n, k; 1), B = [b(i, j)] (0 < i < n'— 1,0 < j < k— 1) represent 
the cyclic DM(i/, k; 1) and C = [c(i, j)] (0 < i < n!, 0 < j < k— 1) represent the cyclic DCA(/c, n' + 1; n'). Recall 
that a(n — 1 ,j) = b{n' — 1 ,j) = c{n',j) = 0 for all 0 < j < k — 1 and a(i, k — 1) = b(i ', k — 1) = c(i", k — 1) = 0 
for all0<i<n — l,0<i' <n' — 1, and 0 < i" < n'. 

Construct a matrix Q = [■ q(i,j)\ where, for 0 < i < n — 1, 0 < i' < n' — 1, 0 < i" < n' and 0 < j < k — 1, 

q(i + i'n,j) = a(i,j) + b(i',j)n, when i ^ n — 1 
q{n — 1 +i"n,j) = c{i" ,j)n. 

Then 

A = {a(i,j) + b(i',j)n — a(i, j') — b(i',j')n \ 0 < i < n — 2,0 < i' < n' — 1} 

U {c(i",j)n — c(i",j')n | 0 < i" < n'}, 

= h n \ {0, n, 2n,..., n(n' — 1)} U {0, n, 2 n ,..., nn'/ 2, nn'/2 ,..., 
n{n' — 1)}. 

Properties PI and P2 follow as in the proof of Lemma T3.51 

□ 


This result can be generalised to construct non-cyclic difference covering arrays as in m- 
We now combine these results with various results of i, m and [20] to obtain results for cyclic DCA(4, 2 p + 
1;2 p), satisfying PI and P2, implying new existence results for row complete 3-MNOLS(2p), where 2 p < 358. 
In the lists below * indicates that the existence of 3-MNOLS(2p) was previously unknown. 

In doing this we will settle the remaining cases in Lemma 13.41 with the exception of 2 p = 146. Here we 
believe that there exists a DCA(4,147; 146) satisfying PI and P2 but have been unable to verify it. 
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Theorem 3 . 9 . {SJ Thm 17.6, p fll] If n is a prime greater than or equal to k, then there exists a cyclic 
DM(n , k ; 1). 

Lemma 3 . 10 . \S(K Lem 2.5] Let n > 5 be prime. Then there exists a cyclic 77DM(4,2n; 2). 

Lemma 3 . 11 . There exists cyclic DCA( 4, 2p + 1; 2 p) for 2p = 50*, 98, 170*, 242, 290, 338. 

Proof. Corollary 13.71 together with Theorem 13.91 and Lemma 13.101 can be used first to construct a cyclic 
HDM(4, 2 p; h) with 2 p(h) = 50(10), 98(14), 170(10), 242(22), 290(10), 338(26) and then the required DCA(4, 2 p+ 
1;2 p). 

□ 


Lemmas 13.13113.15113.171 do not document any new existence results, however they do verify that for the 
given orders cyclic DCA(4, 2p + 1; 2 p) satisfying PI and P2 exist. 

Lemma 3 . 12 . 120 1 Thm 2.1] Let n be an odd positive integer satisfying gcd{n , 9) ^ 3. Then there exists a cyclic 
HDM(4,2n;2). 

Lemma 3 . 13 . There exists cyclic DCA(4, 2p + 1; 2 p) for 2p = 90, 126, 198, 234, 306, 342. 

Proof. Corollary 13.71 together with Theorem 13.91 and Lemma 13.121 can be used first to construct a cyclic 
HDM(4, 2 p; h) with 2 p{h) = 90(10), 126(14), 198(22), 234(26), 306(34), 342(38) and then the required DCA(4, 2p+ 
1 5 2p). 

□ 


Theorem 3 . 14 . 120 . Thm 2.3] Let n > 4 and n = 2“3 /3 p“ 1 .. .p“‘, where {a, /3) ^ (1,0), ai > 0 and the prime 
factors pi >5 for 1 < i < t. Then there exists a cyclic HDM{A , 2n; 2). 

Lemma 3 . 15 . There exists cyclic DCA(A, 2p + 1; 2 p) for 2p = 60, 80, 84, 100, 112, 120, 132, 156, 160, 168, 
176, 180, 204, 208, 224, 228, 240, 252, 264, 272, 276, 300, 304, 312, 320, 336, 352. 


Proof. Corollary 13.71 together with Theorem 13.91 and Theorem 13.141 can be used first to construct a cyclic 
HDM(4, 2p; h) with 2 p{h) = 60(10), 80(10), 84(14), 100(10), 112(14), 120(10), 132(22), 156(26), 160(10), 
168(14), 176(22), 180(10), 204(34), 208(26), 224(14), 228(38), 240(10), 252(14), 264(22), 272(34), 276(46), 
300(10), 304(38), 312(26), 320(10), 336(14), 352(22) and then the required 
DCA(4,2p + l;2p). 


□ 


Theorem 3.16. [20 ': Thm 2-4] Let n = pf 1 .. . p“‘, where cni > 0 and the prime factors pi > 5 for 1 < i < t. 
Then there exists a cyclic HDM(4,4:n;A). 

Lemma 3.17. There exists cyclic DCA(4,2p + 1; 2p) for 2p = 140, 196, 220, 260, 308, 340. 

Proof. Corollary 13.71 together with Theorem 13.91 and Theorem 13.161 can be used first to construct a cyclic 
HDM(4, 2 p; h) with 2 p{h) = 140(28), 196(28), 220(20), 260(20), 308(28), 340(20) and then the required DCA(4, 2 p+ 
1 5 2p). 

□ 


Theorem 3 . 18 . J5J Thm 3.10] A cyclic DM( 3 l , 5; 1) exists for all i > 3. 

Lemma 3 . 19 . There exists cyclic HCt4(4, 2 p + 1; 2 p) for 2p = 216, 270, 324. 

Proof. Corollary 13.71 together with Theorem 13.181 and Theorem 13.141 can be used to first construct a cyclic 
HDM(4, 2 p; h) with 2 p(h) — 216(54), 270(54), 324(54) and then the required DCA(4, 2 p + 1; 2 p). 

□ 

The next result from Yin’s paper |20j is interesting in that it allows us to construct difference covering arrays 
and so nearly orthogonal Latin squares of order 6 n, and gives many values that were previously unresolved (the 
obstruction was the non-existence of MNOLS of order 6). 
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Theorem 3 . 20 . IWl Thru 2.2] Let n be a positive integer of the form p^pf 2 .. .pf*, where cti > 0 and the 
prime factors pi > 5 for 1 < i < t. Then there exists a cyclic HDM(A, 6 n\ 6). 

Corollary 3 . 21 . Let n be an integer of the form p^pf 2 ■ ■ -P^, where on > 0 and the prime factors pi > 5 for 
1 <i <t. Then there exists a cyclic DCA(4,6n + l;6n) satisfying PI and P2. Consequently there exists cyclic 
DCA( 4,2 p + 1; 2 p) for 2 p = 30*, 42*, 66*, 78,102*, 114*, 138*, 150,174*, 

186,210,222,246,258,282,294,318,330,354. 

The following result can be verified using direct constructions given in the later Section |4] 

Lemma 3 . 22 . There exists cyclic DCA(4.,2p + 1; 2 p) for 2 p = 26*, 266, 2 p = 40, 56, 88, 104, 136, 152, 184, 

200, 232, 248, 280, 296, 328, 344 and 2 p = 34*, 58*, 82*, 106*, 130*, 154, 178*, 202, 226, 250, 274, 298, 322, 

346. 

Proof. Corollary 14.51 given in Section U verifies the existence of cyclic DCA(4, 2p + 1; 2 p) with 2 p = 26*, 266. 

Corollary 14.81 given in Section^ verifies the existence of cyclic DCA(4, 2p+ 1; 2 p) with 2 p = 40, 56, 88, 104, 

136, 152, 184, 200, 232, 248, 280, 296, 328, 344. 

Theorem l4.91 given in Section^ verifies the existence of cyclic DCA(4, 2p+l; 2 p) with 2 p = 34*, 58*, 82*, 106*, 130*, 
154, 178*, 202, 226, 250, 274, 298, 322, 346. 

□ 


Lemma 3.23. There exists cyclic DCAfi, 2p + 1; 2 p) for 2p = 24*, 28*, 32, 36*, 44, 48, 52*, 54*. 

Proof. These results have been verified by computer searches. The first column of the DCA(4, 2p + 1,2 p) is 
given by [0,1, 2,..., 2p — 1], the second column by [1,3,..., 2p — 1,2,4,..., 2p — 2] and the third column by 

2 p = 24: [2, 0, 3, 1, 14, 21, 20, 19, 23, 15, 6, 18, 16, 10, 17, 8, 11, 22, 5, 13, 4, 9, 7, 12] 

2 p = 28: [2, 0, 3, 1, 11, 16, 22, 25, 20, 23, 4, 8, 21, 5, 18, 10, 19, 13, 24, 27, 7, 26, 15, 9, 6, 14, 17, 12] 

2 p = 32: [2, 0, 3, 6, 1, 13, 22, 30, 21, 25, 28, 26, 7, 5, 23, 20, 12, 10, 24, 17, 31, 15, 29, 27, 11, 14, 4, 9, 8, 19, 

18, 16] 

2 p = 36: [5, 35, 13, 20, 11, 9, 1, 31, 10, 2, 30, 33, 4, 34, 32, 25, 28, 16, 27, 22, 3, 29, 19, 24, 18, 15, 6, 23, 17, 7, 
0, 8, 14, 12, 21, 26] 

2 p = 44: [39, 13, 26, 21, 35, 3, 17, 16, 40, 28, 38, 25, 6, 10, 34, 5, 18, 30, 43, 15, 19, 36, 7, 24, 32, 14, 4, 0, 31, 
12, 2, 9, 23, 37, 11, 42, 41, 29, 20, 1, 33, 27, 8, 22] 

2 p = 48: [5, 41, 23, 40, 1, 39, 34, 25, 28, 8, 4, 9, 21, 30, 43, 18, 12, 2, 42, 45, 32, 37, 33, 0, 26, 15, 13, 22, 10, 

35, 44, 7, 36, 16, 27, 19, 46, 38, 3, 47, 31, 29, 17, 14, 11, 24, 20, 6] 

2 p = 52: [18, 12, 50, 37, 16, 6, 45, 4, 31, 34, 47, 21, 29, 2, 5, 22, 38, 3, 39, 27, 0, 15, 51, 7, 28, 24, 42, 40, 48, 

32, 9, 26, 20, 11, 1, 41, 19, 35, 43, 13, 49, 33, 14, 17, 46, 8, 36, 23, 10, 30, 25, 44] 

2 p = 54: [6, 5, 31, 27, 20, 38, 19, 4, 30, 51, 3, 52, 49, 14, 48, 23, 41, 12, 25, 0, 32, 40, 21, 50, 9, 45, 16, 1, 46, 
11, 28, 42, 47, 35, 39, 2, 22, 13, 34, 33, 24, 44, 15, 53, 7, 17, 37, 36, 26, 18, 10, 43, 29, 8]. 

□ 


For the remaining values 64, 68, 72, 74, 76, 92, 96, 108, 116, 122, 124, 128, 144, 146, 148, 162, 164, 172, 188, 
192, 194, 212, 218, 236, 244, 256, 268, 284, 288, 292, 314, 316, 332, 348, 356 we were unable to construct cyclic 
difference covering arrays however for completeness and to answer questions about the spectrum we give full 
details verifying existence. It should be noted that it is possible to construct DCA(4, 2p + 1; 2 p) satisfying PI 
and P2 for some of these orders however our construction does not give cyclic difference covering arrays and so 
the details have been omitted here. 

For the Li and van Rees conjecture El Conjecture 5.1] we require two more results from their paper. The 
second result uses group divisible designs: A /C -group divisible design of type g^g^ 2 ■ • •3“'* is a partition Q of a 
finite set V, of cardinality v = Xu=i a i9ii i n t° a i groups of size <?i, 1 < * < s, together with a family of subsets 
( blocks ) B of V such that: 1) if B £ £>, then \B\ £ /C, 2) every pair of distinct elements of V occurs in 1 block 
of B or 1 group of Q but not both, and 3) \Q\ > 1. 

Theorem 3 . 24 . m Thm f.l] Suppose there exists k-MNOLS(2p), k-MOLS(2p), and k-MOLS(n). Then there 
exists k-MNOLS(2pn). 

Theorem 3 . 25 . m Thm 4-5] Suppose there exists a K.-GDD of type g] 1 ... gf 3 . Further suppose that for any 
group size gi there exists a s-MNOLS(gi) and for any block size k £ 1C there exists a s-IMOLS(k). Then there 
are s-MNOLS(Y]i= i <^3*)- 



Lemma 3 . 26 . There exists 3-NMOLS(2p) for 2 p = 76, 92*, 96, 108, 116*, 124*, 128, 144, 148*, 164, 172*, 
188, 192, 212, 236, 244, 256, 268, 284, 288, 292, 316, 332, 348 and 356. 

Proof. The 3-MNOLS(2p), 2 p = 76(12 5 16 1 ), 92(20 4 12), 96(16 5 16 1 ), 108(20 5 8 1 ), 

116(20 5 ,16), 124(20 5 24 1 ), 128(24 5 8 1 ), 144(24 5 24 1 ), 148(24 5 28), 164(28 5 24 1 ), 

172(32 5 12), 188(32 5 28 1 ), 192(32 5 36 1 ), 212(36 5 32 1 ), 236(40 5 36 1 ), 244(40 5 44 1 ), 

256(44 5 36 1 ), 268(44 5 48 1 ), 284(48 5 44 1 ), 288(48 5 48 1 ), 292(48 5 52 1 ), 316(52 s 56 1 ), 

332(56 5 52 1 ), 348(56 5 68 1 ) and 356(60 s 56 1 ) can be constructed applying 5-GDD, that exist by [6j Thm 4.17, p 
258], in Theorem 13.251 Here the bracketed information gives the type of the GDD. □ 

Lemma 3 . 27 . There exists 3-NMOLS(2p) for 2p = 64, 68, 72, 74*, 122, 162, 194, 218 and 314. 

Proof. 3-NMOLS(2p) for 2 p = 64, 68, 72, 74*, 122, 162, 194 and 218 can be constructed by applying 8-GDD(8 s ), 
{7,8,9}-GDD(8 7 6 2 ), 9-GDD(8 9 ), {7,8,9}-GDD(8 7 6 3 ), {7,8}-GDD(16 7 10 1 ), {11,12,13}-GDD(12 12 10 1 8 1 ), 

{11,12,13}-GDD(16 11 10 1 8 1 ), {13,14}-GDD(16 13 10 1 ) and {8,9,10, 11}- 

GDD(32 8 14 2 10 1 ), that exist by finite field constructions, in Theorem l3.251 respectively. □ 

All the results of this section combine to the following theorem. 

Theorem 3 . 28 . There exists a set of 3-MNOLS{2p) for each positive integers p > 3, except possibly p = 73. 
Conjecture 3 . 29 . There exists DGA(4, 2p + 1; 2 p) satisfying PI and P2 of all positive integers p > 3. 

4 Construction of difference covering arrays DCA(4,2m+ 1; 2m) 

This section is devoted to giving new constructions for families of cyclic DCA(4, 2m + 1; 2m) when m = 2k + 1, 
m = 8k + 4 and m = 3k + 2, respectively. In each of these cases the difference covering arrays satisfy PI and 
P2 and so they can be used to construct MNOLS(2m). 

When using a cyclic DCA(4, 2m + 1; 2m) to construct nearly orthogonal Latin squares we strip off the last 
row and last column of zeros. Thus to reduce the complexity of the notation and to avoid confusion, we will 
assume that we are constructing a 2m x 3 array Q = [■ q(i,j )] that satisfies: 

• each column is a permutation of Z-im and 

• A jjf = {q(i,j) — q{i,j')\ | 0 < * < 2to — 1} = {1, 2,..., m, m ,..., 2m — 1}, with repetition retained. 

Also we will use the following notation: q(a, 0) = a (or q(a, 0) = a(a)), q(a , 1) = 6(a) (or g(a, 1) = 6(a)), and 
q{a , 2) = c(a) (or q(a, 2) = c(a)). 

The following lemmas document some well known results, stated without proof, which will be used extensively 
in the proof of subsequent results. 

Lemma 4.1. For all integers x, y, z, gcd(cc + yz, z) = gcd(x , z). 

Lemma 4 . 2 . Let g andp be positive integers and h a non-negative integer. Working modulo 2p, if gcd(g, 2p) = 1 
then 


{gx + h\ 0<x<2p — 1} = Z 2p , 
or if gcd{g , 2 p) = r and h = s mod r then 

{gx + h | 0 < x < 2p/r — 1} = {rx + s | 0 < x < 2p/r — 1}. 

4.1 Construction for general families DCA(4, 2m + 1; 2m) for some odd m 

In this subsection we give a general construction for a difference covering array DCA(4, 2 to + l;2m), for to 
odd. The proof that such a difference covering array exists uses the results presented in the following lemma. 
Note that in this section unless otherwise stated all arithmetic is modulo 2 to. In particular, for i ^ 2 mod 3 a 
non-negative integer, and k = 2i 2 + 7i + 6, we present an infinite family of DCA(4, 2 to + 1; 2to) for to = 2k + 1. 
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Intervals 

h 

12 

13 

Ii 

q(a , 0 ) — a 

[0, m + /] 

[m + / + 1, m — 1J 

[m, rn - f - 1] 

[m — f, 2m — 1J 

q(a, 1) = 6(a) 

a/ + m 

af + m 

(a + 1)/ 

+m — 1 

(a + 1)/ 

+m — 1 

q(a , 2) = c(a) 

_ (a -!)(/ + !) 
-2 

-(a-l)(/ + l) 
+m — 2 

-A/ + i) 

+ra 

~a(f + !) 

6(a) — a 

“(/ - i) 

+m 

“(/ - !) 

+m 

(a + l)(/ —1) 
+ra 

(«+!)(/“ 1) 
+ra 

c(a) — a 

-(a - l)(/ + 2) 
-3 

-(a - l)(/ + 2) 
+m — 3 

—a (/ + 2 ) 
+m 

—a (/ + 2) 

c(a) — 6(a) 

- 0 ( 2 /+ 1 ) 

+/ — 1 + m 

1 

'-h to 

-a(2/+ 1) 

-(/-l) 

-o(2/+ 1) 
-(/ - 1) + m 


Figure 1: Entries are elements of Z 2 TO , where q(a, 0) = a, q(a, 1) = 6 (a) and q(a, 2) = c(a) in the array 
Q = [<?(*> j)]- Rows 5 to 7 give the differences. 


Lemma 4.3. Let f and m be integers such that gcd(f, 2 to) = 2, gcd(f + 2, 2m) = 2, and / 2 + / + 1 = m 
mod 2m. Then 

ged(/,m) = 1, 
gcd(/ + l,m) = 1, 
gcd(/ - l,m) = 1, 
gcd(2/ + 1, m) = 1 

mf = 0 mod 2m 

Proof. Eg 14.11 Note that since / is even, f 2 + f + 1 is odd implying m is odd and hence gcd(/, m) = 1. 
Ea l4.2l Since f+1 = —/ 2 mod m and gcd(/, m) = 1 we have 1 = gcd(/, m) = gcd(/ 2 , m) = gcd(/+l, m). 

EqESl Since / - 1 = (-/ 2 - / - 1) + / 2 + 2/ = /(/ + 2) mod m, gcd(/, m) = 1 and gcd(/ + 2, m) = 1, 

we have 1 = gcd(/(/ + 2),m) = gcd(/ — l,m). 

Eg 14.41 Since 2/ + 1 = —/(/ — 1) mod m, gcd(/,m) = 1 and gcd(/ — 1,m) = 1,1 = gcd(/(/ — 1),m) = 
gcd(2/ + l,m). 

Eo 14.51 This follows from the fact that / is even. □ 

For a suitable choice of /, we divide the domain of a into the subintervals [0, m + /], [m + / + l,m — 1], 
[m, m — / — 1] and [m — /, 2m — 1] where all endpoints are included. 

Example 1. To aid understanding we begin with an example where m = 13 and / = 16 and give the transpose 
of the difference covering array DCA(4, 27; 26). The key to understanding the proof is to recognise that within 
the subintervals I\ = [0,..., 3], I 2 = [4,..., 12], I 3 = [13,..., 22] and I 4 = [23,..., 25], the value of a, 6(a) and 
c(a) increases by a constant “jump”, respectively 1, / = 16 and — (/ + 1) = 9. This implies that the differences 
will also increase by a constant. By carefully choosing the start value on each subinterval it is possible to obtain 
the required values in Z 2 m . The value m = 13 is boldfaced in the differences. 



h 

12 

I 3 

u 


a 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

b( 

a) 

13 

3 

19 

9 

25 

15 

5 

21 

11 

1 

17 

7 

23 

2 

18 

8 

24 

14 

4 

20 

10 

0 

16 

6 

22 

12 

c( 

a) 

15 

24 

7 

16 

12 

21 

4 

13 

22 

5 

14 

23 

6 

0 

9 

18 

1 

10 

19 

2 

11 

20 

3 

25 

8 

17 

6(a) - 

a 

13 

2 

17 

6 

21 

10 

25 

14 

3 

18 

7 

22 

11 

15 

4 

19 

8 

23 

12 

1 

16 

5 

20 

9 

24 

13 

c(a) - 

a 

15 

23 

5 

13 

8 

16 

24 

6 

14 

22 

4 

12 

20 

13 

21 

3 

11 

19 

1 

9 

17 

25 

7 

2 

10 

18 

c(a) — 6( 

a) 

2 

21 

14 

7 

13 

6 

25 

18 

11 

4 

23 

16 

9 

24 

17 

10 

3 

22 

15 

8 

1 

20 

13 

19 

12 
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Theorem 4.4. Let f and m be natui'al numbers such that gcd(/, 2m) = 2, gcd(/ + 2,2m) = 2, / 2 + / + l = m 
mod 2m, and m + 3 < / < 2m — 4. Then a cyclic DCA{A 1 2m + 1; 2m) satisfying PI and P2 exists. 

Proof. The proof is by construction with the values of DCA(4, 2m + l;2m) as given in Figure [Q with the 3 
columns of Q = [q(i,j)\ given by q(a, 0) = a, q(a, 1) = 6(a), q(a , 2) = c(a). We will show that Q has the required 
properties. 

If m + 3 < / < 2m — 4 and / is even then, working modulo 2m, 3 < m + / < m — 4 and so the intervals 
[0,m + /] and [m + / + l,m — 1] are non-empty. Further, m + 4 < m — / < 2m — 3 and so the intervals 
[m, m — / — 1] and [m — /, 2?n — 1] are non-empty. 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 
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Given that / is even and m is odd, by Lemma FT721 


af + m = 1 mod 2 => {b(a) \ a £ I\ U I 2 } = { 2 g +1 \ 0 < g < m — 1 }, 

(a + 1)/ + m — 1 = 0 mod 2 ==> { 6 (a) | a £ I 3 U J 4 } = {2g | 0 < g < m — 1}, 
and so { 6 (a) | 0 < a < 2 m — 1 } = [ 2 m]. 

On each of the subintervals c(a) takes the form ag + h, where g = — (/ + 1), so the “jump” size is —(/ + 1) 
and 

c(m) = 0 , 

c(m + / + 1 ) — c(m — / — 1 ) = —m — /(/ + l) + m — 2 + m 

-/(/+!) —(/ + l)-m 
= -2/ 2 -2/-2-(/ + l) = -(/ + l), 
c( 0 ) - c(m - 1 ) = -(/ + 1 ), 
c(m - /) - c(m + /) = -(/ + 1 ), 
c(m) — c( 2 m — 1 ) = —(/ + 1 ). 

Thus by reordering the subintervals as - 63 , 12 , h, I a, and noting for instance, c(m — / — 1) — (/+1) = c(m + / +1), 
we get {c(a) | 0 < a < 2 m — 1 } = [ 2 m]. 

For 6 (a) — a, c(a) — a and c(a) — 6 (a) we are required to show that for a £ [2m] the differences cover the 
multiset {1,2,..., m — 1, m, m, m + 1,..., 2m — 1} = ([2m] \ {0}) U {m}. 

For 6 (a) — a, the subintervals are taken in natural order / 1 , 12 ,- 63 ,- 64 - Starting at a = 0 and finishing at 
a = 2 m — 1 , we have 6 ( 0 ) — 0 = m = ( 2 m — 1 + 1 )(/ — 1 ) + m = b( 2 m — 1 ) — ( 2 m — 1 ), so the difference m 
occurs twice. Further, the gcd(/ — 1, 2 to) = 1 implies that a(f — 1) + m, 0 < a < m — 1, are all distinct, as are 
(a + 1 )(/ — 1 ) + m, m < a < 2 m — 1 and 

6 (m) — m = (m + 1 )(/ — 1 ) + m = / — 1 , 

6 (m — 1 ) — (m — 1 ) = (m — 1 )(/ — 1 ) + m = —(/ — 1 ). 

Thus there is a jump of —2(/ — 1) between a = m — 1 and a = m and the difference 0 is omitted, implying 
{ 6 (a) — a | 0 < a < 2m — 1} = ([2m] \ {0}) U {m}. 

For c(a) —a, since f + 2 is even and gcd(/ + 2, m) = 1, these values are all distinct on each of the subintervals, 
I /3 U 7i| = m + 1, I /4 U / 2 I = to — 1, and 

c(m) — m = —m(/ + 2 ) + m = m, 
c(m + /) - (m + /) = — (m + / - 1)(/ + 2) - 3 = m, 

(c(0) - 0) - (c(m - / - 1) - (m - / - 1)) = m - /(/ + 2) - 3 = -(/ + 2), 

c( 2 m - 1 ) - ( 2 m - 1 ) = -( 2 m - 1 )(/ + 2 ) = / + 2 , 

c(m + / + 1) - (m + / + 1) = -/ 2 - 2/ + m - 3 = -(/ + 2). 

Thus —(a — 1)(/ + 2) — 3, — a(f + 2) + m = 1 mod 2, hence, 

{c(a) — a | a € /3 U / 1 } = {2g +1 |0<g<m — 1}U {m}. 

In addition, — a(f + 2) + m — 1, —a(f + 2) = 0 mod 2 implies that 

{c(a) — a | a G /4 U I 2 } = { 2 g \ 1 < g < m — 1 }, 

giving {c(a) — a | 0 < a < 2 m — 1 } = ([ 2 m] \ { 0 }) U {m}. 

For c(a) — 6 (a), since gcd(2/ + 1, 2m) = 1, these values are all distinct on the subintervals, c(m + / + 1) — 
6 (m + / + 1 ) = — 2/ 2 — 2 / — 2 + to = m = m + 2/ 2 + 2 / + 2 = c(m — / — 1 ) — 6 (m — / — 1 ), and 

c(0) — 6(0) — (c(m — 1) — 6(m — 1)) = —(2/ + 1), 

c(m + /) - 6(m + /) = 2/ + 1, 

c(m - /) - 6(m - /) = -(2/ + 1). 

Thus when the subintervals are reordered to I 2 , h, Ia, h we may verify that {c(a) — 6 (a) | 0 < a < 2 to — 1 } = 
([2m] \ {0}) U {to}. Note that the values of c(a ) — 6 (a) start and finish on m and the value 0 is omitted between 
a = to + / and a = m — f. □ 
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Corollary 4.5. Let i ^ 2 mod 3 be a non-negative integer, k = 2i 2 + 7i + 6 and m = 2k + 1. Then there exists 
an infinite family of cyclic DCA(4, 2 m + 1; 2to) ’s satisfying PI and P2. 

Proof. Taking / = m + 3 + 2i, then / is even. In addition to + 3 < / < 2m — 4, since i > 0 and 2 m — 4 = 

(to + 3) + (to — 7) = (to + 3) + (4 i 2 + 14* + 6) > to + 3 + 2i = f. Now 

f 2 + / + 1 = 4* 2 + 14* + 13 mod 2m 

= 2(2* 2 + 7* + 6) + 1 = to mod 2m. 

Further, applying Lemma l4~Tl repeatedly. 

gcd(/, 2 to) = 2(gcd(2* 2 + 8* + 8, (2* 2 + 8* + 8) + 2* 2 + 6* + 5)) 

= 2(gcd(2* + 3, 2i 2 + 4i + 2 + (2* + 3))) 

= 2(gcd(2* + 3, 2(* + l) 2 )) = 2(gcd(2* + 3, * + 1)) = 2 


Also, 


gcd(/ + 2,2m) = 2(gcd(2* 2 + 8* + 9, (2 i 2 + 8* + 9) + 2 i 2 +6i + 4)) 

= 2(gcd(2* + 5, 2(* + 2)(* + 1))) 

= 2(gcd(2* + 5, (* + 2)(* + 1))). 

Now gcd(2* + 5, * + 2) = gcd(2(* + 2) + 1, * + 2) = 1. Whereas 

gcd(2* + 5, * + 1) = gcd(2(* + 1) + 3, * + 1)) = gcd(3, * + 1) 

1 when * + 1 = 0 mod 3 or equivalently i = 2 mod 3. 

Thus taking * ^ 2 mod 3 we can construct a DCA(4, 2 to + 1; 2to) as per the Theorem 14.41 □ 

4.2 Construction of difference covering arrays DCA(4,4?ri + l;4m) 

In this subsection we give a general construction for a difference covering array DCA(4, 4 to + 1; 4to). It will be 
shown that for all non-negative integers k , such that 3f (2 k + 1), this construction gives an infinite family of 
DCA(4,16fc + 9; 16fc + 8). The proof that such a difference covering array exists uses the results presented in 
the following lemma. Note that in this section unless otherwise stated all arithmetic is modulo 4 to. 

Lemma 4.6. Let f and to be natural numbers such that to = 2 mod 4, gcd(/, 


4m) = 2, gcd(/ — 1, 4 to) = 1, and f 2 + f — 2 = 2m mod 4m. Then 

gcd(2m + 2 — /, 4 to) = 4, 

(4.6) 

gcd(2m — / + 1, 4 to) = 1, 

(4.7) 

gcd(2?7* — 2/ + 2, 4 to) = 2, 

(4.8) 

to/ = 2to mod 4 to. 

(4.9) 


Proof. Eg 14.61 Rewriting 2 to + 2 — f = f 2 + f — 2 + 2 — f = f 2 mod 4 to and assuming gcd(/, 4 to) = 2 gives 

gcd(/ 2 , 4 to) = 4. 

EofTTl Since 2 to — / + 1 is odd, the gcd(2?T* — / + 1,4m) is odd. Assume there exists an odd x such that 

x|4to and x\(2 to — / + 1), then x\m and so x \(/ — 1). But the gcd(/ — 1,4 to) = 1, so x = 1. 

Ea 14.81 Assume that there exists x such that x\{m — f + 1) and x\2m. Since to — / + 1 is odd, x is odd 
and so x\m. Consequently x\ (/ — 1) and a;|4m, implying x = 1. 

Ea 14.91 It follows that to/ = to(2to — / 2 + 2) = 2m 2 — to/ 2 + 2to = 2to mod 4m. □ 

Theorem 4.7. Let f be a natural number and m = 4 k + 2, where k is a non-negative integer, such that 

gcd(/, 4 to) = 2, gcd(/ — 1,4m) = 1, and f 2 + / — 2 s 2 to mod 4 to. Then a cyclic DCA[4,,Am + l;4m) 
satisfying PI and P2 exists. 
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Intervals 

h 

12 

I 3 

Ii 

q(a , 0) = a 

[0, m — lj 

[m, 2m — 1] 

[2m, 3m — 1] 

[3m, 4 m — 1] 

q(a , 1) = 6(a) 

(a + l)/-l 

(a + l)/-l 

af 

af 

q(a, 2) = c(a) 

(a+l)(2m-/ + 2) 
-1 

a(2m — / + 2) 

—m 

(a+l)(2m-/ + 2) 
+m — 1 

a(2m — / + 2) 

6(a) — a 

(«+!)(/-!) 

(0 + 1 H/- 1 ) 

«(/ - !) 

a (/ - i) 

c(a) — a 

(a + l)(2m-/ + l) 

a(2m — / + 1) 
—m 

(a+l)(2m-/+l) 

+m 

a(2m — f + 1) 

c(a) — 6(a) 

(a + l)(2m - 2/ + 2) 

a(2m - 2/ + 2) 
— / — m + 1 

a(2m - 2/ + 2) 
+3m - / + 1 

o(2m - 2/ + 2) 


Figure 2: Entries are elements of Z 4 m , where q(a, 0) = a, q(a, 1) = 6(a) and q(a, 2) = c(a ) in the array 

Q = WJ)}- 


Proof. The proof is by construction with the values of DCA(4,4?n + l;4m) as given in Figure [2 with the 3 
columns of Q = [q(i,j)\ given by q(a, 0) = a, q(a, 1) = 6(a), q(a , 2) = c(a). We will show that Q has the required 
properties. 

For 6(a), since / is even, Lemma 14.21 implies that 

(a + 1)/ — 1 = 1 mod 2 =>• {6(a) | a € I\ U 7 2 } = {2 g + 1 | 0 < g < 2m — 1}, 
af = 0 mod 2 ==> (6(a) j a 6 4 U 74} = {2g 1 0 < g < 2m — 1}, 


and {6(a) | 0 < a < 4m — 1} = Z 4 m . 

For c(a), since gcd(2m — / + 2,4m) = 4, Lemma IT21 implies that 


- / + 2) + m 

- 1 

= i 

mod 

4 = 

=> {c(a) | a 

G M = 

= {4ff 

+ 1| 

0 

< 






g < m — 

1}, 





a(2m — / + 2) 

= 0 

mod 

4 = 

=> {c(a) | a 

G /4} = 

= {4ff 

|o< 

g 

< 






m — 1}, 






2m - / + 2) 

- 1 

= 3 

mod 

4 = 

=4- {c(a) | a 

G7!} = 

= {4ff 

+ 3| 

0 

< 






g<m- 

1}, 





!m — / + 2) - 

- m 

= 2 

mod 

4 = 

^ {c(a) | a 

G 7 2 } = 

= {4ff 

+ 2| 

0 

< 






g < m — 

!}• 






Thus the set of values {c(a) | a € Z 4 m } = Z 4 m . 

For 6(a) — a, c(a) — a and c(a) — 6(a) we are required to show that for a € Z 4 m the differences cover the 
multiset {1,2,..., 2m — 1, 2m, 2?n, 2m + 1,..., 4m — 1} = \ {0}) U {2m}. 

For 6(a) — a, the gcd(/ — 1,4m) = 1, and 


6(2?n) — 2 m = 2m(f — 1) = 2m, 

6(2m — 1) — (2m — 1) = (2m — 1 + 1)/ — 1 — (2m — 1) = (2 m)(/ — 1) = 2m, 
6(4m — 1) — (4m — 1) = —(/ — 1), 

6 ( 0 ) — 0 = / — 1 . 


So using a “jump” of / — 1 and ordering the subintervals as I 3 , 14 , I\, I 2 we obtain the difference 2m twice and 
the difference 0 is omitted between a = 4m — 1 and a = 0 implying that {6(a) — a | 0 < a < 4m — 1} = 
(Z 4m \ {0}) U {2m}. 

For c(a) — a, the gcd(2?n — / — 1,4m) = 1, and 

c(m) — m = m(2m — / + 1) — m = 2m, 
c(3m — 1) — (3m — 1) = (3m — 1 + l)(2m — / + 1) + m 
= 2m, 

c(3m) — 3m — (c(2m — 1) — (2m — 1)) = (m + l)(2m — / + 1) + m 

= 2m — / + 1, 
c(0) — 0 = 2m — / + 1, 

c(4?n — 1) — (4m — 1) = (4m — l)(2m + / — 1) 

= -(2m-/ + 1), 

c(2m) — 2m — ( c(m — 1) — (m — 1)) = (m + l)(2m — / + 1) + m 

= 2m — / + 1. 
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So using a “jump” of 2m — / + 1 and ordering the subintervals as 4 2 ,4 4 , we obtain the difference 2 m 
twice and the difference 0 is omitted between a = 4m — 1 and a = 0, implying \c(a) — a I 0 < a < 4m — 1} = 
(Z 4 m \ {0}) U {2m}. 

For c(a ) — b(a), and 

c(3m) — 6(3m) = m(2m — 2/ + 2) = 2m, 
c(m — 1) — b(m — 1) = (m — 1 + l)(2m — 2 + 2) = 2m, 
c(4m — 1) — 6(4m — 1) = — ( 2m — 2/ + 2), 
c(0) - 6(0) = 2m- 2/+ 2, 
c(2m) — 6(2m) — (c(2m — 1) — 6(2m — 1)) = 2m — 2/ + 2. 

Then since 

2m — 2/ — 2 = 0 mod 2 ==> {c(a) — 6(a) | a € / 4 U / 1 } = {2 g | 1 < g < 

2m — 1} U {2m}, 

—/ + 1 = 1 mod 2 => {c(a) — 6(a) | a G / 2 U 73 } = {2<? + 1 | 0 < g < 

2m — 1 }, 

implying {c(a) — 6(a) | 0 < a < 4m — 1} = (Z 4m \ {0}) U {2m}. □ 

Corollary 4.8. For k > 0 such that k ^ 1 mod 3 a cyclic DCA(4,4m + l;4m) satisfying PI and P2 can be 
constructed as described in Theorem 

Proof. Given m = 4fc + 2, take / = 2m — 2. Then / = 8 k + 2, and 

gcd(/, 4m) = 2(gcd(4/c + 1, 8A: + 4)) = 2(gcd(4/c + 1, 2(4fc + 1) + 2)) = 2. 

In addition 

gcd(/ — 1,4m) = gcd(2?n — 3,4m) = gcd(8fc + 1,16fc + 8) 

= gcd(8/c + 1, 2 k + 1) since 2{(8k+l) 

= gcd(6A:, 2k+ 1) = 1, if 3{(2k + 1). 

Also 


f 2 + / - 2 = 64fc 2 + 32fc + 4 + 8 k + 2 - 2 = Ak{l 6 k + 8) + 8 k + 4 = 2m mod 4m. 

Hence, / = 2m —2 satisfies the assumptions of Theorem 077] and we can construct a DCA(4,16/c + 9; 16k+ 8) 
for k such that 3{(2fc + 1) as described in this theorem. □ 

4.3 Construction of difference covering arrays DCA(4,2?n + l; 2m), where m = 3^ + 2 

In this subsection we give a general construction for a difference covering array DCA(4,2m + l;2m), where 
m = 3/r + 2. The proof that such a difference covering array exists uses the result presented in the following 
lemma. Note that in this section unless otherwise stated all arithmetic is modulo 12fc + 10, k > 0. 

Theorem 4.9. Let yi be an odd positive integer. Then there exists a cyclic HC14(4,6/x + 5; 6/i + 4) satisfying 
PI and P2. 

Proof. Since p > 1, n = 6 p + 4 > 10, since 3 is prime, gcd(3, 6 p + 4) = 1. Hence gcd(3/r + 4, 6 p + 4) = 
gcd(3/x,6^ + 4) = 1. 

Let p = 2k + 1, k > 0, then 

p( 8 p + 2) = (2k + l)(6fc + 5) = 3/i + 2, and 
(3p + 2) 2 = 3p + 2. 

That is, (n/2 ) 2 = n/2 mod n. 
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Intervals for 

h 

I 2 

h 

a 

[0, M — 1] 

[/4, 2/i] 

[2/4 + 1, 3/4 + 1] 

q(a , 0) = a(a) 

3a -|- 3/i + 4 

3a+ 2 

3a 4" 3/i 4” 4 

</(a, 1) = 6(a) 

3a(/i —|— 1) —|— 2/i -|- 2 

3a(/i 4~ 1) 4~ 2/i -|- 2 

3a(/i 4~ 1) 4" 2/i 4- 2 

</(a,2) = c(a) 

a(3/i —|— 4) —|— 5/i + 4 

a(3/i H - 4) 4~ 5/i 4" 4 

a(3 (i 4" 4) 4" 5/i 4~ 4 

6(a) — a(a) 

3a/4 — /4 — 2 

3a/4 + 2/4 

3a/4 — /4 — 2 

c(a) — a(a) 

a(3/i + 1) + 2/i 

a(3/i 4" 1) 4" 5/i 4" 2 

a(3/i 4-1) 4" 2/i 

c(a) — 6(a) 

a -|- 3/i + 2 

a 4" 3/i 4" 2 

a 4" 3/i 4" 2 

Intervals for 

h 

h 

is 

a 

[3/i + 2,4/i + 2] 

[4/i + 3, 5/i 4~ 2] 

[5/4 + 3,6/4 + 3] 

S' 

e 

II 

o' 

3a -|- 3/i 4- 3 

3a 4~ 1 

3a 4" 3/i 4~ 3 

?(a, 1) = 6(a) 

3 a(/i + 1) + 2/i + 1 

3a(/i 4-1) 4- 2/i + 1 

3a(/i 4" 1) 4" 2/i 4-1 

II 

CnT 

a(3/i H - 4) -|- 5/i + 4 

a(3/i 4" 4) 4" 5/i H - 4 

a(3/i 4" 4) 4" 5/i 4~ 4 

6(a) — a(a) 

3a/i — /4 — 2 

3a/r + 2/4 

3 a/4 — /4 — 2 

c(a) — a(a) 

a(3/i + 1) + 2/i 1 

a(3/i 4" 1) 4" 5/i 4- 3 

a(3/i 4~ l) 4" 2/i 4" 1 

c(a) — 6(a) 

a -|- 3/i + 3 

a 4" 3/i 4" 3 

a 4~ 3/i 4" 3 


Figure 3: Entries are elements of 2^+4, where q(a, 0) = a(a), q(a, 1) = 6(a) and q(a, 2) = c(a) in the array 

Q = [<?(*,j)]- 


The proof is by construction with the values for DCA(4, 6 fj, + 5; 6/4 + 4) as given in Figure G2 with the three 
columns of Q = [q(i, j)] given by q(a, 0) = a(a), q(a, 1) = 6(a) and q(a , 2) = c(a). 

Since gcd(3, n) = 1, {3a |0<a<n — 1} = Z„, by Lemma, 14.21 Further a( 3/4 + 2) = 3(3/4 + 2) + 3/4 + 3 = 1 
and a(2/j,) = 6/4 + 2 and there is a “jump” of 3 between a(4/r + 2) and a(0); a(/4 — 1) and a(5/4 + 3); a(6/4 + 3) 
and a(2/i + 1); a(3/4 + 1) and a(4/r + 3); a(5/4 + 2) and a(/r), respectively. Thus reordering the subintervals as 
gives (a(a) |0<a<n-l} = Z„. 

For 6(a), 

3a (/4 + 1) + 2/4 + 2 = 0 mod 2 ==+ { 6 (a) | « € Ji U 62 U 13 } 

= {25 I 0 < g < 3/4 + 2} 

3 a (/4 + 1) + 2/4 + 1 = 1 mod 2 =>• { 6 (a) | a g 64 U J 5 U Ig} 

= + 1|0<3^3/x+1} 

For c(a), since gcd( 3/4 + 4, 6/4 + 4) = 1 Lemma l4~2l implies that {c(a) | 0 < a < 6/4 + 3} = Z 6 /i+ 4 . 

For 6(a) — a(a), since gcd(3/4, 6/4 + 4) = 1, 

b(n) - a(p.) = 3/4 + 2, 

6(4/4 + 2) — u(4/4 + 2) = 3/4 T 2, 

6(5/4 + 3) — a(5/4 + 3) — (6(2 /r) — a( 2/r)) = 3/4 
6(4/4 + 3) — a(4/4 + 3) — (6(/i — 1) — a(/4 — 1)) = 6/i 
6(2/i + 1) — a(2/i + 1) — (6(5/t + 2) — a(5/i + 2)) = 3/4 
6(0) — a(0) — (6(6/i + 3) — a(6/i + 3) = 3/4. 

So using a “jump” of 3/t and ordering the subintervals as I 2 , I&, Ii, I 5 , i 3 , I 4 , we obtain the difference 3/4 + 2 
twice and since there is 6/4 between 6 (a) — a(a) for a = /i — 1 and a = 4/4 + 3 the difference 0 is omitted implying 
that { 6 (a) — a(a) | 0 < a < 6/4 + 3} = (Z„ \ {0}) U {n/2}. 

For c(a) — 6 (a), since gcd( 3/4 + 1, 6/4 + 4) = 2, we have 
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and 


c(3/i H- 2) — ci(3fi -|- 2) 

c(m) - a(/z) 

c(0) - a(0) 
c(M - 1) - a(n - 1) 


c(5/i + 3) — a(5/t + 3) 
c(2/t) — a(2/t) 
(c(6/t + 3) — a(6/t + 3)) 
(c(4/t + 2) u(4/4 + 2)) 

c(2/i + 1) — a(2/i + 1) 
c(5/i + 2) — c(5/i + 2) 
(c(3/t + 1) - a(3/i + 1)) 
(c(4/t + 3) - a(4// + 3)) 


3/r + 2, 

3/i + 2, 

— (3/i + 3) 
-(3/i + 3) 
3/i + 1 
3/i + 3 
-(3/i + 3) 
-(3/i + 3). 


Reordering the intervals as Iq , 1 . 4, /2 and 13, fi, Is and using a regular “jump” of — (3/i + 3), with the jump of 
6 /i + 2 between a = 2/i + 1 and a = 5/i + 2, being the exception, we have the difference 3/t + 2 twice and the 
difference 0 omitted, thus |c(a) — a(a) | 0 < a < 6/1 + 3} = (Z 6ai+ 4 \ {0}) U {3/i + 2} with repetition retained. 
For c(a) — 6(a), we note that 


c(0) -6(0) = 3/i + 2, 
c(3/i + 1) - 6(3/i + 1) = -1, 

and so the values of c(a) — a(a) on the subinterval Ii U I 2 U I 3 cover the set {3/i + 2,..., —1}. Also 


c(3/i + 2) — 6(3/i + 2) = 1, 
c( 6/1 + 3) — 6 ( 6/1 + 3) = 3/i + 2, 

and so the values of 6(a) — c(a) on the subinterval I 4 U I 5 U Iq cover the set {3/i + 2, ..., 6/1 — 1}. Consequently 
(6(a) — c(a) | 0 < a < 6/1 + 3} = ([ 6/1 + 4] \ {0}) U {3/i + 2} with repetition retained. □ 


4.4 Infinite families 

The construction of Theorem [OJ] constructs sets of three MNOLS of orders 10, 22, 34,46 mod 48. The construc¬ 
tion of Corollary 14.81 constructs sets of three MNOLS of orders 8,40 mod 48. Combined with the constructions 
of [7], there is a construction of three MNOLS for 8,10,14,22, 34, 38,40,46 mod 48. There are infinite families 
constructed from Corollary 14.51 and from results of Li and van Rees [IT], but these cannot be described mod 48. 
It is an open question as to why 48 features in many of the constructions. 
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